Gates for the Kane Quantum Computer in the Presence of Dephasing 
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In this paper we investigate the effect of dephasing on proposed quantum gates for the solid-state 
Kane quantum computing architecture. Using a simple model of the decoherence, we find that the 
typical error in a CNOT gate is 8.3 x 10 -5 . We also compute the fidelities of Z, X, Swap, and 
Controlled Z operations under a variety of dephasing rates. We show that these numerical results 
are comparable with the error threshold required for fault tolerant quantum computation. 



I. INTRODUCTION 

One of the most exciting advances in physics has been 
the development of quantum algorithms 0,0 which out- 
perform their best known classical counterparts. These 
algorithms are described in the absence of noise, and de- 
coherence. In experiment this will certainly not be the 
case. In this paper we investigate numerically how a sim- 
ple model of decoherence affect gates on the Kane quan- 
tum computer Q . The Kane quantum computer is one of 
a number of promising silicon based quantum computer 
proposals 0, |H IE 0, Si 

Dephasing in systems similar to the Kane architecture 
has been investigated since the introduction of the spin 
echo technique [9j. The nuclear spin and electronic spin 
decoherence times of P donors in Si is relatively long [T(J, 
[H [H E Q E E H3 ■ For example, recently a time 
of 60ms was measured for the electronic dephasing time, 
E3- Although dephasing times are comparatively 
long, if left unchecked the accumulated errors introduced 
by dephasing will destroy coherence in the computation. 

Using quantum error correction protocols it may be 
possible to correct the errors caused by decoherence. 
El El El El E! To successfully reduce the overall 
error in the system, we must correct errors faster than 
they accumulate, in a fault-tolerant manner 1241 . This 
consideration leads to an error threshold [25l I26tl27l E3 • 
Typically such a threshold requires the probability of in- 
troducing an error in each gate to be below 1 x 1CU 4 
to as low as 1 x 10~ 6 . For the Kane architecture the 
exact threshold is still under investigation [2^. In this 
paper we ask, if it is possible for the Kane architecture 
to achieve this error threshold. To do this we must know 
how much error is introduced by each of the gates used 
in the Kane quantum computing architecture. 

This paper shows the results of simulations of gates 
on the Kane architecture in the presence of dephasing. 
Single qubit gates presented here are similar to those 



used in Nuclear Magnetic Resonance (NMR) |3jj, |3jJ for 
rotations of individual qubits. Voltage fluctuations on 
the 'A' gate and stochastic modeling of the system has 
been investigated analytically in Refs. [32L l33l |. The 
two qubitgates presented here use non-adiabatic pulse 
schemes [34j . 

The analysis given here is a direct analogue to that of 
Fowler et al js3 f° r the adiabatic CNOT gate. The gates 
analyzed in this paper are simpler, faster, and potentially 
higher fidelity than the adiabatic gate j34|. These gates 
do not rely on complicated pulse shapes, but simply turn- 
ing on or off the voltages applied. In contrast to adiabatic 
gates, the timing of these gates could easily be run off a 
digital clock cycle. In addition, whereas the adiabatic 
CNOT is required to be applied up to three times to 
create an arbitrary two qubit gate, using non-adiabatic 
schemes it is possible to create an arbitrary two qubit 
gate directly. For example, the swap gate analyzed here 
would require three adiabatic CNOT gates to construct. 
Using non-adiabatic schemes we are able to construct it 
in a single pulse sequence |34j . which is much faster, sim- 
pler and higher fidelity than the corresponding adiabatic 
scheme. 

We compare each of the gates analyzed to the error 
threshold for fault tolerant quantum computing. We sim- 
ulate the master equation for typical values of spin de- 
phasing expected in the Kane architecture. We find that 
the error in the gates analyzed is less than or comparable 
to that required for fault tolerant quantum computation. 

This paper is organized as follows. In Section [H] we 
describe the simple model of decoherence which we use 
and present the master equation for the system. Sec- 
tion lllll presents the results for one qubit gates, including 
free evolution in subsection IIII Al Z rotations in subsec- 
tion llll Bl and X rotations in subsection llll CI Two qubit 
non-adiabatic gates are shown in Section llVI These in- 
clude the CNOT gate in subsection IIV Al and the swap 
gate and controlled Z gates in subsection IIV Bl Finally, 
conclusions are drawn in Section IV! 
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II. THE MASTER EQUATION 

A brief introduction to the Kane quantum computing 
architecture [j| is given here. The Kane architecture con- 
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sists of P donor atoms embedded in Si. The orientation 
of the nuclear spin of each P donor represents one qubit. 

When placed in a magnetic field applied in the z di- 
rection, Zeeman splitting occurs. This is given by the 
Hamiltonian 

Hb = -gnV-nBZn + fi B BZ e , (1) 

where Z is the Pauli Z matrix, and the subscripts e (n) 
indicate electronic (nuclear) spin. The magnetic field, 
B, may be controlled externally. The application of a 
resonant rotating magnetic field adds the following terms 
to the spin Hamiltonian: 

H ac = -gnVnBae [X n COs(uj ac t) + Y n Sm(u ac t)} 

+^BB ac [X e cos(oj ac t) + Y e sin(u; ac t)] . (2) 

The electronic spins couple to their corresponding nuclear 
spins via the hyperfme interaction, 

H A = Aa e ■ cr n , (3) 

where A is the strength of the interaction. The design 
of the Kane quantum computer calls for control of the 
strength of the hyperfme interaction externally by ap- 
plying appropriate voltages to 'A' gates. The electronic 
spins couple to adjacent electrons via the exchange inter- 
action 

Hj = Jcr ei ■ cre 2 , (4) 

where e\ and ei are two adjacent electrons, and J is the 
strength of the exchange interaction which may be con- 
trolled externally through the application of voltages to 
the 'J' gates. 

Altogether the spin Hamiltonian of a two donor system 
is given by 

2 

H s = H Bi + H Ai + Hj + H aCi . (5) 

i=l 

The times and fidelities of the gates naturally depend 
on exactly which parameters are used to calculate them. 
For many of the gates in this paper the typical parameters 
shown in Table^were used. These parameters are similar 
to the parameters used for the pure state calculations in 
0. 

A simple model of decoherence was used for these cal- 
culations. There are many different decoherence mech- 
anisms, but our model only considers pure dephasing 
(without cngergy relaxation) . Whereas dephasing is cer- 
tainly not the only source of decoherence, it likely to 
be the dominant effect on a time scale shorter than the 
energy relaxation (dissipation) time, T±. For example, 
Feher and Gere [lj] measured T\ n > 10 hours for nu- 
clear spin at a temperature of T = 1.25K, B = 3.2T and 
T\ e f» 30 hours under similar conditions. In contrast, 
experimentally measured times for Ti have been much 
shorter. Gordon and Bowers measured T 2e = 520/is 



Description 


Term 


Value 


Unperturbed Hyperfme 
Interaction 


A 


0.1211 x 10~ 3 meV 


Hyperfme Interaction 
During Z Rotation 


A z 


0.0606 x 10~ 3 meV 


Hyperfme Interaction 
during X Rotation 




0.0606 x 10" 3 meV 


Constant Magnetic Field 
Strength 


B 


2.000T 


Rotating Magnetic Field 
Strength 


B ac 


0.0025T 


Hyperfine Interaction 
during Interaction 


Au 


0.1197 x 10" 3 meV 


Exchange Interaction 
during Interaction 


Ju 


0.0423 meV 



TABLE I: Typical parameters used for numerical calculations. 



for P : Si at T = 1.4K in isotopically enriched 28 Si. Chiba 
and Harai have also measured the electronic decoher- 
ence times of P : Si, finding a rate of = lQOfis. For 
the nucleus, recent results for the nuclear spin of a 29 Si 
nucleus show a maximum value of T"2„ = 25s j3|| . 

Recently Tyryshkin et al [l8[ obtained an experimen- 
tal measurement of T^^ = 14.2ms at T = SAK and 
T2 c = 62ms at T = 6.9K for a donor concentration of 
0.87 x 10 15 cm -3 in isotopically pure Si. At millikelvin 
temperatures the decoherence time is likely to be even 
longer. Additionally these measurements were carried 
out in a bulk doped sample, and in our case we will 
be considering a specifically engineered sample. Inter- 
actions such as exchange and dipole-dipole interactions 
contribute to the coupling between electrons. In some ex- 
periments, such as in Ref. these potentially benefi- 
cial coupling have been treated as sources of decoherence, 
but in the operation of a quantum computer these inter- 
actions can be either be decoupled or used to generate 
entanglement, useful for quantum computation. Hence, 
it is expected that T2 c may even be longer than those re- 
ported in 0] . Nevertheless, we use the value of 60ms as a 
conservative estimate for electronic dephasing time. We 
expect the nuclear dephasing times to be several orders 
of magnitude bigger than electronic dephasing times. We 
choose the following parameters to be typical of P : Si the 
systems we are considering: 

T 2c = 60ms, (6) 
T 2n = Is. (7) 

The typical errors presented in the tables contained in 
the next three sections are evaluated at these typical de- 
phasing times. 

The simple decoherence model we consider corresponds 
to the master equation 

p = ~[H s ,p]-C[p}, (8) 
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where the dephasing terms are given by 



In the rotating frame, the master equation has the so- 
lution 



C[p] 



T e [Z ei ,[Z ei ,p]\ +T n [Z ni ,[Z ni ,p\] (9) 



Characteristic dephasing rates, and T2 n are related 
to the dephasing rates by the equations: 



T 2c 
T 2n 



1 
1 

4Tv 



(10) 
(11) 



We define fidelity (and therefore error) in terms of the 
actual state after applying an operation, p, and the in- 
tended state after that operation, p' . Due to systematic 
errors and decoherence these states will not necessarily 
be the same. When comparing against a pure state //, 
the fidelity F of an operation is defined as 



F(p,p')=Tv(pp'); 
Error is defined in terms of fidelity 

E(p,p / ) = 1-F(p,p'). 



(12) 



(13) 



Typically we would like to know the greatest error pos- 
sible for any input state. This is a computationally dif- 
ficult problem. In the results which follow, the approach 
taken is to calculate the fidelity for each of the compu- 
tational basis states, and each of the input states which 
would ideally generate a Bell state. This has two main 
benefits. The first is that a high fidelity indicates that 
the gate is successfully creating or preserving entangle- 
ment. The second is that Bell states are superposition 
states, which are susceptible to dephasing. We also calcu- 
lated the effect of each gate on the four Bell input states 
for the CNOT gate. For typical parameters, Bell input 
states give similar fidelities to those shown in this paper. 

Throughout this paper we will use the states |0) and 
|1) to represent the nuclear spin up and spin down states 
respectively. We will use the | |) and | |) to represent 
electronic spin up and spin down states respectively. 



III. ONE QUBIT GATES 

A. Free Evolution 

The spin of an isolated nucleus undergoing Larmor pre- 
cession in the presence of a magnetic field, B [37| • In this 
case it is easy to solve the master equation with dephas- 
ing exactly. Considering only the nuclear spin, we have 

H s =g n fi n BZ n . (14) 
The decoherence terms has only the single term 

c[ P ] = r n [z n ,[z n ,p\\, (15) 

= 2T N (p- Z nP Z n ). (16) 



P(t) = 



Poo(0) poiWe- 4 ^ 4 
Pw(0)e- 4rNt pn(0) 



(17) 



This has the effect of exponentially decaying the off di- 
agonal terms of the density matrix, but leaves the diago- 
nal components unchanged. For a single isolated nuclear 
spin, the simple model has no effect on eigenstates of 
Z n (i.e., there is no relaxation process for these states). 
In contrast, it has a dramatic influence on superposition 
states whose off diagonal terms decay exponentially (i.e., 
dephasing). Two such states are 



1 

71 

1 

71 



(|o>- 

(10) 



We can easily calculate the expectation value of the Pauli 
X matrix, (X), for the |+) state 



(X) = Tr(A», 

= exp(-4rAr<). 



(18) 
(19) 



For a single nuclear spin coupled to an electronic spin 
via the hyperfmc interaction the Hamiltonian is given by 



H S =H B + H A . 



(20) 



We assume that electron is initially polarized by the large 
magnetic field, B. The evolution of this Hamiltonian 
was calculated for their typical values [given in Eq. 10, 
Eq. (ZJ and Table |J. The fidelity after different times 
is shown in Fig. ^ This figure shows the Bloch sphere 
radius, given by 



r= ^(X)2 + (Y)2 + (Z)2. 



(21) 



A pure state has a radius of one, and a radius of less than 
one indicates a mixed state. An initial state of | J, +) 
was used. The radius decays at a rate governed by the 
nuclear decoherence time, which in this case is Ti n = Is. 
As is expected, this decay is the same as the well known 
solution to the Bloch equations. 



B. Z Rotations 

Z rotations on the Kane architecture may be performed 
by varying the Larmor precession frequency of a single 
qubit |34j . Graphs showing the error in the Z gate at 
different dephasing rates are shown in Figs. Eland EH 

Figure [21 shows the error at different rates of decoher- 
ence for both electrons, T e , and nuclei, r„, for a single 
qubit in the | J, 0) state. The calculated error does not 
significantly depend on the electronic or nuclear dephas- 
ing rate. The pure | J, 0) state is not affected by decoher- 
ence terms in the master equation. Therefore the only 
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FIG. 1: Effect of dephasing on the free evolution of a single 
embedded P atom in the Kane architecture. 



10" 




r N (■-) 



FIG. 2: Error in the Z gate for | J. 0} initial state at differing 
rates of dephasing. 



State 


Systematic Error 


Typical Error 


|0> 

l+> 


3.8 x l(T b 

1.9 x 10 -6 


3.8 x l(T b 

1.9 x 10 -6 


Maximum 


3.8 x 10 _e 


3.8 x l(T b 



TABLE II: Summary of Z gate error. 



effect of dephasing occurs when the hyperfine interaction 
rotates this state. The effect of dephasing on this state 
is negligible. 

The error in the Z operation for an initial state of | j 0) , 
is primarily due to systematic error of 3.8 x 10~ 6 . This 
error is due to the hyperfine interaction coupling between 
electrons and nuclei. This allows a small probability of 
finding the electrons in an excited state. At typical rates 
of decoherence for the short duration of a Z rotation (ap- 
proximately 21ns), systematic error is the dominant ef- 
fect. For the | j 0) state, the typical error is 3.8 x 1CP 6 . 

In the previous section we noted that superposition 
states, | J. +) and | J, — ), arc affected by dephasing terms 
more than eigenstates of Z . This is illustrated in Fig. 
Electronic dephasing times have little effect on the 
overall fidelity. As the nuclear dephasing rate increases, 
the fidelity decreases, with a maximum error of approx- 
imately 0.5. This indicates all quantum coherence was 
lost and we are in classical mixture of the states | J, 0) 
and | I 1). For typical rates of dephasing [given in Eq. 
©, and Eq. J7J], the error is found to be 1.9 x 1CT 6 
which is largely due to systematic error. 

The maximum error, for typical rates of dephasing, 
of any of the states tested for the Z gate is 3.8 x 10~ 6 . 
This error is largely due to systematic effects rather than 
dephasing. This error suggests it is theoretically possible 
to do a Z rotation with an error of less than the 1 x 10~ 4 
limit suggested for fault tolerant quantum computing. 
The results for the Z gate are summarized in Table ITT1 
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FIG. 3: Contour plot showing the error in the Z Gate for 
| J, +) initial state evolving at differing rates of dephasing. 
Contour lines connect, and are labelled by, points of equal 
error. 



C. X Rotations 

X rotations are performed using a resonant magnetic 
field, B ac . The error in the X rotation was found at dif- 
ferent rates of dephasing. This is shown for the two basis 
states in Figs. 0] and El Evident on both of these graphs 
is a valley which levels out at a minimum error. This er- 
ror is primarily due to systematic error in the gate. For 
the | J, 0) initial state the systematic error is found to 
be 2.3 x 1CP 6 , and for the | J, 1) initial state the sys- 
tematic error is found to be 4.9 x 10~ 6 . As the nuclear 
dephasing rates, T2 n , increases, the fidelity of the opera- 
tion decreases. The fidelity of the operation drops to 0.0 
indicating that, in the limit of large dephasing rates, the 
rotating magnetic field does not have the desired effect 
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State 


Systematic Error 


Typical Error 


10) 
|1> 


2.3 x l(T b 
4.9 x 1(T 6 


3.8 x l(T b 
6.4 x 10 -6 


Maximum 


4.9 x l(T b 


6.4 x 10"° 



TABLE III: Summary of X gate error. 
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FIG. 4: Contour plot showing the error in the X gate for | J. 0) 
initial state evolving at differing rates of dephasing. Contour 
lines connect, and are labelled by, points of equal error. 



of a resonant magnetic field. 

In Figs. 0] and [SI we see that the fidelity of the X ro- 
tation depends weakly on the electronic dephasing rate, 
r 2e . The solutions to the equations under these condi- 
tions show that the principal cause of error is the electron 
becoming excited to a higher energy level. 

Under typical conditions, the | J. 0) initial state has 
an error of 3.8 x 10~ 6 . The | J, 1) state has an error of 
6.4 x 10~ 6 . In this operation dephasing has a much more 
important role than in Z rotations. One reason for this 
is that an X rotation takes longer, approximately 6.4/is. 

The error in the X gate is summarized in Table II I II 
The maximum error from the two basis states tested for 
the X gate was 6.4 x 10~ 6 . The error induced in this 
operation is less than a threshold of 1 x required for 
fault tolerant quantum computing. 



IV. TWO QUBIT GATES 
A. The CNOT Gate 



The CNOT gate is specified by 



T 1 X = 



10 
10 
1 
10 



(22) 
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FIG. 5: Contour plot showing the error in the X Gate for 
j I 1} initial state at differing rates of dephasing. Contour 
lines connect, and are labelled by, points of equal error. 



which may be created using the steps specified in Ref. 
[3^ . In the following discussion of two qubit gates, if 
not explicitly stated, all initial electron spin states are 
assumed to be | J. J,). 

The error in the CNOT in the presence of dephasing 
was found by numerically solving the master equation 
[Eq s. © and |JS}] for the appropriate pulse sequence 
[34j . For different rates of dephasing, different fidelities 
are obtained. Fidelities were calculated for each of the 
four computational basis states, and for the evolution 
leading to the four Bell states. Fidelities for the four 
computational basis states are shown in Fig. |SJ 

Each of the states shows a minimum error when the 
rates of electronic and nuclear dephasing arc slow. The 
remaining error is due to systematic error in the gate. 
Some sources of systematic error for the CNOT oper- 
ation include off-resonant effects, excitation into higher 
electronic energy levels, and imperfections in the pulse 
sequences (such as the break-down of the second or- 
der approximations used to derive appropriate times for 
pulse sequence HU)- Systematic error for each of the 
states |00), |01), |10) and |11) are 4.0 x 10" 5 , 2.6 x 10~ 5 , 
1.9 x 10~ 5 , and 2.9 x 10 -5 respectively. For evolution 
starting in an initial Bell state, we find the systematc er- 
ror is 3.5 x 10~ 5 , 3.4 x 10~ 5 , 1.9 x 10" 5 and 2.6 x 10~ 5 . 
Sytematic error for states resulting in a Bell state are 
shown in Table llVl 

In each of the four states, as the dephasing rate in- 
creases, the fidelity decreases. In the limit of large de- 
phasing rates, the computational basis states tend to stay 
in their original states. This is particularly evident from 
the graphs of 1 00) and |01) which have higher fidelity 
(lower error) at high dephasing rates. In contrast, the 
states 1 10) and |11) have lower fidelities (highter error) 
at high rates of dephasing. For example, the 1 00) state 
stays in the 1 00) state after the CNOT operation. At such 
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(c) 1 10) initial state 

FIG. 6: Contour plots showing error in the CNOT operation 
lines connect, and are labelled by, points of equal error. 




(d) 1 11) initial state 



different rates of electronic and nuclear dephasing. Contour 



high rates of dephasing, not even the single qubit rota- 
tions described in Section IIII CI apply, which form part 
of the CNOT gate operation. For these states quantum 
coherence has been lost. When we consider, for example, 
the state \tp) = 1/V2(|00) + |01)) we find that at a high 
dephasing rate of r„ = T e = 54 x 10 6 s _1 that the error 
of the gate is 0.5. In this case, quantum coherence has 
been lost between the two states, and the qubit evolves 
to a completely mixed state. 

Electronic dephasing rates play a much bigger role in 
two qubit gates than in single qubit gates. 

For the typical dephasing times, T2 e and T2 n [given in 
Eq. ©, Eq. Q], we find the error for the states |00), 
|01), 1 10) and |11) are 8.3 x 10" 5 , 6.8 x 10~ 5 , 6.2 x 10~ 5 , 
and 7.2 x 10 -5 respectively. The error for an initial state 



of one of the four Bell states is found to be 6.0 x 10 5 , 
6.0 x 10" 5 , 4.4 x 10" 5 , and 5.1 x 10" 5 . Typical errors 
for states resulting in a Bell state are shown in Table 
IIVI This implies that under our very simple decohcrence 
model, the maximum error in the CNOT gate in any basis 
state is 8.3 x 10~ 5 . This is only marginally under the 
threshold of 1 x 10~ 4 required for fault tolerant quantum 
computation. 

Errors for each of the computational basis states are 
shown in Table llVl and the maximum error of any of the 
four computational basis states or states with evolution 
leading to a Bell state is plotted in Fig. [7| 

These results are directly analogous to calculations 
made for the adiabatic CNOT gate [3f|. We have used 
the same noise model as was used in their numerical sim- 
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State 


Systematic 


Error 


Typical Error 


|00> 


4.0 x 10 




8.3 x 10" 


b 


|01) 


2.6 x 10 




6.8 x 10" 


5 


|10) 


1.9 x 10 


~ B 


6.2 x 10" 




111) 


2.9 x 10 




7.2 x 10" 


5 


|00) + 11) 


2.9 x 10 




7.0 x 10" 


5 


|00) - |11) 


3.2 x 10 


— 5 


7.3 x 10" 


5 


|01) + |10) 


3.1 x 10 


— 5 


7.2 x 10" 


5 


|01) - |10) 


2.3 x 10 


— 5 


6.4 x 10" 


5 


Maximum 


4.0 x 10 


— b 


8.3 x 10" 


5 



TABLE IV: Summary of CNOT gate fidelities. 
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FIG. 7: Contour plot of the maximum error in basis and 
Bell output states of the CNOT gate shown as a function of 
electronic and nuclear dephasing rates. Contour lines connect, 
and are labelled by, points of equal error. 



ulations. In calculations for the adiabatic CNOT gate at 
dccohcrence times, for electronic and nuclear dephasing 
rates of T2 e = 200/xs and = 10s respectively, the 
maximum error of any of the four basis states for the 
adiabatic CNOT gate was found to be 'just over 10~ 3 ' 
|35j . In comparison, for the same conditions, we find that 
the maximum error in the non-adiabatic gate is 7 x 10~ 3 . 
Under these conditions both non-adiabatic and adiabatic 
gates give similar fidelities. 

Where does this error come from at the dephasing 
rate specified above, for the nonadiabatic CNOT gate? 
The error in an X rotation under these conditions is 
3 x 10~ 6 and the error in the entangling operation of 
the gate U m (j), is 4 x 10~ 3 and therefore it is clear 
that the two qubit entangling operation is the major 
source of error. When electronic dccohcrence times are 
short, any electron mediated operations will be affected 
by this decohcrence. By increasing the exchange interac- 
tion strength, the time required for the electron mediated 
operation may be reduced. For example, at a strength of 
J = 0.0529mcV the error decreases to 1 x 1CP 3 and the 
maximum error in the CNOT gate is also 1 x 10~ 3 for 
any of the computational basis states. 



The advantage of nonadiabatic gates over adiabatic 
gates is that the pulse schemes required are much sim pler , 
faster, and as at the conditions considered in Ref. (35| . 
the two schemes have approximately the same fidelity. 
Considering the electronic decoherence times measured 
in Ref. 0| it is likely that electronic dephasing rates are 
not as large as considered in Ref. [3f|. At lower rates 
of dephasing, we approach the systematic error, which 
may be smaller for nonadiabatic gates than adiabatic 
gates [331 . Another distinct advantage of the nonadia- 
batic gates is that any two qubit gate may be made. This 
allows us to construct two qubit gates directly (such as 
the swap gate), which are faster and higher fidelity than 
expressing them as combinations of CNOT gates and sin- 
gle qubit rotations. If the Kane computer was being run 
from a digital clock cycle, non-adiabatic two qubit gates 
could be controlled at discrete times, and do not require 
the continuous and sophisticated pulse shapes required 
for adiabatic gates operating at this speed and fidelity. 



B. The Swap Gate and Controlled Z Gate 

Similar calculations to those calculated for the CNOT 
gate were carried out for the swap gate and the controlled 
Z gate. The swap gate is specified by the matrix 

10 
10 
10 
1 



(23) 



The controlled Z gate is specified by the matrix 



T 1 Z = 







-1 



(24) 



The circuit which may be used to create the swap gate 
may be found in Ref. |34j . 

The master equation was solved numerically for each 
pulse sequence. The error in the swap gate was calculated 
for each basis state, and each state whose output state is a 
Bell state. Note that for these two gates the states which 
give Bell states as output are themselves Bell states. A 
separate simulation was completed for each combination 
of nuclear and electronic dephasing times. The maximum 
error of any of the basis states has been plotted in Fig. 
[H] for the swap gate, and Fig. [5] for the control Z gate. 

Similar features that were evident for the CNOT gate 
are visible in these figures. The corresponding errors are 
shown in Table IVI for the swap gate, and Table IVI1 for the 
control Z gate. 



V. CONCLUSION 

In conclusion, we have investigated the effect of de- 
phasing on the Kane quantum computing architecture. 
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State 


Systematic 


Error 


Typical Error 


|00> 


3.9 x 10 




9.0 x 10" 




|01) 


1.4 x 10 




7.9 x 10" 


-5 


|10) 


1.6 x 10 


~ 5 


8.0 x 10" 




111) 


3.8 x 10 




8.9 x 10" 


_5 


|00) + 11) 


5.3 x 10 




1.4 x 10" 


-4 


|00) - 111) 


7.4 x 10 


— 5 


1.6 x 10" 


-4 


|01) + |10) 


1.7 x 10 


— 5 


1.0 x 10" 


-4 


|01) - |10) 


1.5 x 10 


— 5 


1.0 x 10" 


-4 


Maximum 


7.4 x 10 


— b 


1.6 x 10" 


-4 



TABLE V: Summary of swap gate error. 




FIG. 8: Contour plot of the maximum error in basis states 
and Bell output states of the swap gate shown as a function of 
electronic and nuclear dephasing rates. Contour lines connect, 
and are labelled by, points of equal error. 



We used a simple model of decoherence and investigated 
how this model affected proposed gate schemes on the 
Kane quantum computer. For typical decoherence rates 
[given in Eq. © , Eq. {7|)] , these results are summarized 
in Table IVTT1 

Each of the errors, for typical rates of dephasing, found 
here are close to the error threshold required for fault 
tolerant quantum computation. If the temperature is 
lowered, and coupling between qubits is not considered 



State 


Systematic 


Error 


Typical Error 


|00) 


3.0 x 10 


— b 


8.1 x 10" 


ta 


|01) 


1.5 x 10 


— 5 


6.6 x 10" 


5 


|10) 


1.1 x 10 


— 5 


6.2 x 10" 


5 


111) 


3.8 x 10 


-5 


8.9 x 10" 


4 


|00) + |11) 


3.6 x 10 


-b 


9.4 x 10" 


5 


|00) - 111) 


3.3 x 10 


— 5 


9.2 x 10" 


5 


|01) + |10) 


1.7 x 10 


— 5 


7.5 x 10" 


5 


|oi> - |io> 


1.1 x 10 


— 5 


7.0 x 10" 


5 


Maximum 


3.8 x 10 


— 5 


9.4 x lO- 


5 



TABLE VI: Summary of controlled Z gate error. 




FIG. 9: Contour plot of the maximum error in basis states 
and Bell output states of the Y\Z gate shown as a function of 
electronic and nuclear dephasing rates. Contour lines connect, 
and are labelled by, points of equal error. 



Gate 


Typical Error 


Systematic 


Error 


Time 


Z 


3.8 x 10~ B 


3.8 x 10 


-« 


0.02/is 


X 


6.4 x 10" 6 


4.9 x 10 


-6 


6.4/xs 


CNOT 


8.3 x 10~ 5 


4.0 x 10 


-5 


16.0/is 


Swap 


1.6 x 10" 4 


7.4 x 10 


-5 


19.2/xs 




9.4 x 10~ 5 


3.8 x 10 


-5 


16.1/is 



TABLE VII: Times and errors for each of the gates investi- 
gated. 



a decoherence process as in Ref. [18|, it is likely that 
the typical decoherence times for the Kane architecture 
may be further reduced, and therefore unambiguously 
under the threshold required for fault tolerant quantum 
computation. 

Construction and operation of the Kane quantum com- 
puter is extremely challenging. In the actual physical 
system there will undoubtedly be noise and decoherence 
processes not considered in our simple physical model. 
This substantial effort would never be able to achieve 
its ultimate goal of a working quantum computer if there 
were fundamental reasons why such a computer could not 
operate. In this paper we investigated the one such effect 
on proposed gates for the Kane quantum computer. Our 
simulations indicate that errors due to dephasing, the 
dominant form of decoherence in the Kane architecture, 
do not rule out fault tolerant quantum computation. 
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